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Abstract 

Recent  results  have  shown  that  several  H 2  and  H2- 
related  problems  can  be  formulated  as  convex  pro¬ 
grams  with  a  finite  number  of  variables.  We  present 
an  interior  point  algorithm  for  the  solution  of  these 
convex  programs  and  illustrate  its  application  with 
the  standard  LQR  design. 

1.  Introduction 

It  has  been  shown  recently  that  a  number  of  H2  and 
fU-related  problems  can  be  formulated  as  convex  pro¬ 
grams  with  a  finite  number  of  variables  —  quadratic 
stabilization  [1],  mixed  B’ljBca  and  multicriterion 
LQG  problems  (see  [2]  and  references  therein).  The 
common  idea  underlying  these  results  is  that  though 
the  original  problem  is  not  convex,  a  clever  change  of 
variables  [3]  makes  it  convex. 

In  this  paper,  we  present  a  systematic  procedure  for 
transforming  the  convex  programs  resulting  from  H2- 
related  problems  above  into  optimization  over  Affine 
Matrix  Inequalities.  We  then  present  a  simple  inte¬ 
rior  point  method  for  their  solution.  Though  our  pre¬ 
sentation  is  through  the  simple  LQR  design  example, 
the  techniques  readily  extend  to  the  more  complicated 
problems  cited  above. 

2.  The  LQR  problem 

Consider  the  linear  time  invariant  system  described 
by  the  state  equations 

x  =  Ax  +  Bu  +  w 


'  Ri  0 

u 

.  0  Q1-  _ 

X 

where  u  is  the  control  input,  w  is  unit  intensity  white 
noise  and  z  is  the  output  signal  of  interest.  The  LQR 
problem  is  to  design  a  feedback  controller  from  the 
state  x  to  the  control  input  u  which  minimizes  the 
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H 2  norm  between  w  and  z  [4].  It  is  known  that  the 
optimal  feedback  law  is  a  constant  state  feedback  u  = 
—Kx  and  optimizes  the  following  program: 

min  Tr (QP)  +  Tr (R1/2KPKT  R1/2) 
subject  to 

(A  -  BK)P  +  P(A  -  BKf  + 1  <  0  and  P  =  PT  >  0. 

By  defining  a  new  quantity  Y  =  KP ,  the  above  prob¬ 
lem  can  be  written  as 

mm  Tr(QP)  +  Tr(R1/2YP-1YTR1/2)  (2) 

subject  to 

AP  +  PAt  -BY  -  YtBt  +  I  <  0  and  P  =  PT  >  0, 
which  is  a  convex  program  (see,  for  example,  [1],  [2]). 


3.  Transformation  to  Optimization  over 
AMI’s 

The  general  structure  of  the  convex  program  (2)  is 

min  J(P,Y) 

(p,Y)eA 

where  A  is  some  convex  constraint  set,  and  J  a  perfor¬ 
mance  index.  We  show  how  this  can  be  transformed 
into  the  problem 

min  7 

C(7,Z)>  0 

where  (7,  Z)  is  a  new  set  of  variables  and  C  is  sym¬ 
metric  and  an  affine  matrix  function  of  (7  ,Z).  The 
inequality  C'(z)  >  0  is  called  an  Affine  Matrix  In¬ 
equality  (AMI). 

Example:  The  LQR  problem 

The  objective  function  of  program  (2)  consists  of 
the  sum  of  two  terms.  It  is  easily  shown  that  the 
second  term 

$(P,  Y)  =  T y{RiI2YP~1Yt R112) 
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can  be  expressed  as 

<S>(P,Y)  =  min(Tr(Y)) 

'  X  R^2Y  1 

YT R}!2  P  >  ' 


Then  let 


Ci(7,  P,  Y,  X )  :=  -Tr(QP)  -  Tr(X)  +  7, 


C2( 7,  P,  Y,  X)  :=  -AP  -  PAt  +  BY  +  YT BT  -  I, 

v  /?! /2\y 

C3(7,P,Y,X):=  [  yTRl/2  p 
C'( 7,  P,  Y,  X)  :=  diag(CT,  C2,  C3). 


The  optimization  problem  (2)  can  now  be  written 


:=  Y*, 

Compute  a  lower  bound  7 ; . 

end 

Remark  1:  An  initial  lower  bound  7 ;  can  be  chosen 
to  be  0. 

Remark  2:  Computing  the  analytic  center  of  a 
convex  bounded  set  (a_center(diag(C,  70  —  7)  >  0)) 
needs  an  initial  point  interior  to  the  constraint. 
(yu,  Zu)  is  such  a  point. 

Remark  3:  The  lower  bound  7 ;  is  computed  from 
7*  and  the  Hessian  of  the  barrier  function  of  the  con¬ 
straint  expressed  at  (7 *,  Y*).  For  more  information, 
see  [5],  [6]  and  also  [7]  (these  proceedings). 


minimize  7 
C(1,P,Y,X)>0 


(4) 


5.  Conclusion 


which  indeed  is  of  the  form  (3). 

Well-Posedness 

We  say  that  the  convex  program  (3)  is  well-posed  if 
for  every  real  7  the  set  {Z  \  C'( 7,  Z )  >  0}  is  compact 
or  empty.  We  observe  the  following  without  proof: 

Proposition  1:  The  program  (4)  corresponding  to 
the  LQR  problem  is  well-posed  if  R  is  positive  defi¬ 
nite,  (A,  B)  is  controllable  and  ( Q,A )  is  observable. 

Under  similar  assumptions,  all  the  other  problems 
cited  in  the  introduction  enjoy  the  same  property. 

4.  Computational  aspects 

Problem  (3)  is  a  convex  non-differentiable  opti¬ 
mization  program,  and  the  ellipsoid  algorithm  or  Kel¬ 
ley’s  cutting-plane  algorithm  [4]  may  be  used  to  solve 
it.  Recently,  the  work  of  Nemirovski  et  al.  has  led  to 
the  development  of  interior  point  algorithms  based  on 
the  notion  of  the  analytic  center  for  a  set  of  convex 
constraints  [5];  these  algorithms  seem  to  hold  great 
promise. 

We  will  describe  one  such  interior  point  algorithm, 
called  the  method  of  centers.  Given  an  initial  feasible 
point  (yu,Zu)  for  constraint  C  in  program  (3),  and 
a  desired  absolute  accuracy  e  on  the  optimum,  the 
algorithm  is  as  follows: 

while  yu  -  7;  >  e, 

7o  :=  7 u  +  e, 

(7*,  Z*)  :=  a_center  (diag(C,  70  -  7)  >  0), 

In  ■=  7*  , 


Through  the  LQR  example,  we  have  outlined  a  sys¬ 
tematic  procedure  for  transforming  convex  optimiza¬ 
tion  programs  arising  from  ^-related  problems  into 
optimization  over  AMI’s.  We  have  also  briefly  de¬ 
scribed  a  simple  interior  point  method  for  their  solu¬ 
tion.  Our  procedure  easily  applies  to  problems  in  [1] 
and  [2],  and  more  generally  to  many  quadratic  Lya¬ 
punov  function  shaping  problems. 
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